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ABSTRACT 

This paper describes an extended-image tracking technique based on the maximum likelihood estimation. The target 
image is assume to have a known profile covering more than one element of a focal plane detector anay. ^assumed 
that the relative position between the imager and the target is changing with tune and the received target image has 
each of its pixels disturbed by an independent additive white Gaussian noise. When a rotation-mv^t i^ment 
between imager and target is considered, the maximum likelihood based image tracking technique described in this 
paper is a closed-loop structure capable of providing iterative update of the movement estimate by calculating e 
toop feedback signals from a weighted correlation between the currently received target mage and the previously 
estimated reference image in the transform domain. The movement estimate is then used to duect the imager ^o 
closely follow the moving target. This image tracking technique has many potential apphcations, including free-space 
optical communications and astronomy where accurate and stabilized optical pointing is essen 1 . 

Keywords: Optical pointing; Extended- image tracking; Image stabilization; Maximum likelihood estimation; Opti- 
cal communications 


1. INTRODUCTION 

Free-space laser communication is a promising technology to meet the demand of cost-effective broadband support for 
future P NASA missions by allowing spacecraft communication subsystems with lower mass and power consumption 
In order to prepare for development and validation of various deep-space laser communication technologies. JPL 
begun construction at its Table Mountain Facility in California of a 1-meter telescope ‘oLe 

Communication Telescope Laboratory and capable of tracking spacecraft from low Earth orbit to d^p space _One 
of these vital technologies is the acquisition and tracking of laser signals to ensure accurate porting of the receiving 
telescope to the laser transmitter onboard the spacecraft for an optical space-to-ground link (referred to as an optical 
downlink) which may be extended as long as several astronomy units. For those deep-space apphcations, it is not 
uncommon that the required pointing accuracy is on the order of microradians to ensure the qualify ^ of mi opfcad 
downlink in the presence of many factors that cause pointing error. For example uncompensated P latf ° r “ 
and/or jitter, sensor noise and/or bias, and atmospheric propagation effects, such as image dancing, blurring, and 
can all contribute to the error in optical pointing. 2 ’ 

This paper presents an optical downlink acquisition and tracking technique for the toe-pointmg "“J” 1 
based optical terminal. Figure 1 shows the functional block diagram of a typical receiving system of a ground-based 
o^al SStaTft includTa wrangle acquisition-aid telescope which is aligned with a much larger mam e k^pe 
uLd to receive the downlink laser signals. The maligned telescopes are driven by a coarse-pomting controUer winch 
is capable of keeping the targeted spacecraft within the FOV of the main telescope. The downlink laser beam coming 
off the mirror assemblies of the main telescope is then directed to the optical detector by a two-axis f ^ 8t ^ 
mirror. By controlling this fast-steering mirror through a separate fine-pointing controller, the incoming 
can be tracked and the laser spot is centered at the optical detector to ensure optimal reception. 

. Th e research described in this paper was carried out by the Jet Propulsion Laboratory, California Institute of Technology under 
contract with the National Aeronautic and Space Administration. 




Figure 1. The ground receiving system. 


The commonly used fine-pointing techniques involve the laser image acquisition and tracking on the focal plane 
array and can be categorized into two types: the centroid-type algorithm and the correlation-type algorithm. The 
former calculates the image centroid solely from the detected laser intensity profile, however, the latter performs a 
detailed feature comparison between the detected intensity profile and the reference profile which may be pre-stored 
or synthesized. In general, the centroid-type techniques work best if the acquired laser spot is small and has high a 
signal-to-noise ratio. On laser spot extended several pixels on the focal plane array, the correlation-type techniques 
perform better, even with low signal-to-noise ratios. The proposed technique is based on a correlation-type extended- 
source images acquisition and tracking technique developed for optical transceivers onboard the spacecraft. To 
adopt the extended-source model for the downlink when the laser beam is usually characterized as from a point- 
source is justified by the fact that, due to atmospheric propagation, diffiraction and other system imperfections, the 
laser spot on the focal plane can cover more than one pixel of the detector array like an image of an extended-source. 
Furthermore, for data detection, it is a common practice to intentionally de-focus the incoming laser beam to ensure 
that at least a portion of the detector array can capture the signal in the presence of optical turbulence. 

The development of this extended-image acquisition and tracking scheme is based on the maximum likelihood es- 
timation, where the uncertainties between the received image and the reference image are modeled as independent 
additive white Gaussian disturbances. It is also assumed that the relative position between the spacecraft and the 
ground-based optical terminal is changing with time, but only in a rotation-invariant fashion on a 2-D plane. As 
suggested by this study, the received laser image obtained from a detector array is first correlated in the trans- 
form Hnnnitin with the reference image derived from the priorly known or synthesized intensity profile of a distant 
point source. The coordinate of the onboard optical transmitter is then estimated and tracked using the developed 
open-loop acquisition and closed-loop tracking algorithms, respectively. The optimal acquisition requires solving two 
nonlinear equations to estimate the coordinate. A suboptimal estimate exists by solving linearized version of the 
maximum likelihood criterion when computation complexity becomes a concern. A closed-loop tracking algorithm 
motivated by the maximum likelihood criterion is developed as well for continuously tracking of the translation 
movement between the spacecraft and ground terminal. The closed-loop tracking structure has its loop feedback 
signals formulated as weighted transform-domain correlations between the received laser image and the previously 






sub-pixel resolutions in a high disturbance environment. 

In this paper .notion 2 provide. . geper.1 deecription of the mrthemrttad model Md the dtacrete flamer tmnrform 
of mi image 'uS the tedeograpUc repremntation. The effect of rotatkm-invmtmit tr.mdrt.on 

research efforts in section 5. 


2. MATHEMATICAL BACKGROUND 
2.1. Representation of an Extended Image 

The image of incoming laser beam detected by an M x N focal plane array at time t„ *“ * 

represented by a sum of the source image, s,(m,n), and the random disturbance, n,(m,n), as follows 

r i (m,n) = ai(m,n) + n J (m,n), m = 0, - 1, n = 0,1 ,. .., N - 1. (1) 

With an additive white Gaussian random disturbance model, n,(m,n) is assumed to be an independent zenvmean 
Gaussian random variable with variance of for all m and n. 

The discrete Fourier transform of the received image at time U becomes 

■R.i{m,n) = Si(m,n) +M(m,n) 

where the transform-domain source image and random disturbance are 


( 2 ) 


and 


p=0 q=0 

M(m,n) = M £ N £n l (p, q )e-M*’+*<) 


(3) 


(4) 


2.2. Translation Movement 

related to the previous image at U by /g\ 

9i+i(m,n) = si(m — x i, n - Vi) ' ' 


In the transform domain, it can be easily shown that the relationship becomes 

5j + i(m, n) = Si(m,n)e 

where 


„ / m n \ 
e m ,n,i = -2ir \jj x l + JjVl ) 


( 6 ) 


(7) 




image. 



3. MAXIMUM LIKELIHOOD ESTIMATE OF TRANSLATION VECTOR 


Base on the Gaussian assumption stated in Eq. (1), the maximum likelihood estimator will declare the estimated 
translation vector if 

p(rt + 1 | xi,y t ;si) = max p(f}+i | xj, p,; si) (8) 

where, for notational convenience, rj+i and s) are the vector representation (known as the lexicographic form) of the 
corresponding received image matrix r; + i(m, n) and the reference image matrix 5j(m,n), respectively and 


P(ri + i I Xj,yj;s)) = 


(V2n<7i) 


1 — t^tIKi»i -fi‘i II J 


(9) 


is the conditional probability density function of fj+i given that the translation vector is (xj, j/i) and the reference 
at the beginning of this movement is sj. The maximum likelihood criterion stated in Eq. (8) is equivalent to 


min || ri+i - Jj +1 || 3 = min || fj+i - L„, *,{$)} || a =|| f J+ i - L fl>fil {«i} || 2 


( 10 ) 


where L z , y {*} is defined as a translation operator which moves the operand by a translation vector (x,y). It has 
been shown 4 that the maximum likelihood criterion stated in Eq. (10) can finally be reduced to 


max 



|u; i (m,n)|cos(^ mint / - 6 , 



=11 n+i - Lfj.fi, {3} 


(11) 


where 

wi(m,n) = ft|+i(m,n)$*(m,n) = |u>J(m,n)|e <^m ’“• , (12) 

is the pixel-by-pixel product of the transform-domain received image 7£*+i(m,n) and the complex conjugate of the 
transform-domain reference image Si(m y n ). Note that the likelihood function to be maximized in Eq. (11) can be 
rewritten as 

{ M—l N-i >| 

jjjj 53 ^2'R'i+i(m,n) s ;(m,n)e- ie — 1 \ (13) 

m=0 n— 0 ) 

where Re{-} represents the real part of a complex quantity. It is clearly indicated that the likelihood function involves 
the average over all pixels of the pixel-wise multiplied received and reference images in the transform domain, as well 
as the phase to be estimated. 

Taking the partial derivatives of Eq. (13) with respect to xj and yi and equating them to zero, we have 


M-lN-l 

53 53 m|u;i(m,n)|sm(£ m , niJ - 0 m> „,i ) = 0 

m«0 n**0 


M-lN-l 

53 53 n l ty K m > n )l sin (£m, n,» - &m,n,l ) = 0 
m=0 n=0 


(14) 


which form a set of nonlinear equations to be solved for the maximum likelihood estimates of xj and yj. 

When the extended-source is close to being acquired, the phase differences, (£ m ,r»,J - 0m,n,j), are small and the 
approximation of sin(x) « x can be applied to Eq. (14), rendering the suboptimal linear estimator 


M-liV-1 M- 1 N-l 

53 53m|w,(m,n)|U.n. J = 53 53 ™M m , n )|0m.n,J 

m =0 n =»0 m *0 n = 0 

M-1AT-1 M-l N~l 

EE n|u; t (m,n)|$ m> „,j = 53 53 n M m ’ n >i*"* ,n,< 

m *0 n »0 m **0 n *0 


( 15 ) 




Figure 2. The Open-loop Maximum-Likelihood Extended Image Acquisition. 

The estimated vector (ft, ft) obtained from solving this suboptinml liaem mmdmum likelihood criterion most tatiaiy 

-S E E m|u>i(m 1 n)|( mA i = ft £ E ">(*■"» + »EE 


m=0 n=0 

M-1N-1 M-1N-1 

2tT -H. f“i m=0 n=0 m=0 "=° 


m=0 o=0 

M-1N-1 


m=0 n=0 

M-l Af-1 


(16) 


m*=0 n=0 ~ ' 

lit a distant star. Also note that the transform-domain images can be realized by usmg a single ens. 


4. IMAGE TRACKING LOOP 

domain received image K l+ i (m,n) and the translated reference image 

£i+i( m > n ) — <S|(m, n)e** m ■"• , 

from previoos estimstioo to be continuously monitored. The pixel-ertta Product of R, t i(m.n) and *?„<">.») 
be expressed by 

Cl+i{tti, n) = 72e| + i(m,n)5*(m,n)e ** m,f *’ 1 

= |S((m,n)| a e*"*"‘ + M+i (m,n)5f (m,n)e-*”*"‘ ( 18 ) 

where 6 m , n ,i is the estimate of 6 m , n ,i in Eq. (7) and the estimation error is 

= -2tt g(*. - ii) + jj{yi - y«)] 

- - 2 * (S A * + £ A *) 

where A* and A v are the associated estimation errors. 


( 19 ) 






According to the mAvlmnm likelihood criterion derived in the previous section, the real part of Eq. (18) shall be 
maximized over the entire detector array, rendering the likelihood function 

53 H Re {Cj+i (m,n)} (20) 

m«0 n~0 

which is the correlation between ft i+1 (m,n) and £j+i (m,n). It turns out without surprise that £i and yi can be 
obtained by solving a set of two nonlinear equations formed by setting the partial derivatives of Eq. (20) with respect 
to xi and yi to be zero. However, in deriving the tracking algorithm to continuously update the estimates, two loop 
feedback signals are similarly formed as the partial derivatives of Eq. (20) with respect to A z and A v , rendering 
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- M-1AT-1 M-l N—l 
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m*0 n**0 


m*0 n«0 


( 21 ) 


(22) 


where 

^/m.n) = ^-Re{M + i(m,n)5r(m,n)e- rf —'} 

= ^-Rc{M+i(m,n)«S,*(m,n)e-" — '} 

become the effective noises in the loop operation. Equations (21) and (22) characterize the relationship between 
the estimate errors, A x and A y , and the loop feedback signals, £ x £y However, to solve for A x and A y from 
these nonlinear equations can be a quite challenging task. With a reasonable assumption valid when the phase error 
remains small during the tracking mode, one can substitute Eq. (19) for sin(^ m , n ,j) in Eqs. (21) and (22). The 
resulting simultaneous equations are linear for A x and A v and can be easily solved, yielding 


A (7nE[£ x ] — C mn K[€y 

x ~ C m C n -Cl 


Ay = 


mn 

Cmn E^] - C m E[s y ] 
*2 

mn 


~ CmQn 


(23) 

(24) 


where £[•] denotes the statistical expectation and 
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are coefficients that can be calculated from the transform-domain reference image of the previous iteration at tj. 




Figure 3. The Extended-Image Tracking Loop. 

An extended-image tracking loop structure can be realized based upon the above discussion and is depicted in 
pfg t Th^tranlrm-Sn received image {**!<«, »)} * first multiplied pixel-wise wxth a prepay transU^ 
trLform-domain reference image established according to the estimate fo, yi) from the previous iknkn< ^ 
r g averaged over the entire detector array, the correlation result is used to compute the loop faadtek «j 
and s v . It ta important to note that, in real implementation, these two loop feedback signals are calculated from 

Cj+i(m,n) as the following weighted sums 


£x = 1 E mIm{C i+ i(m,n)} 
M m=0 n=0 
M-1N-1 

£„ " S 51 «Im{Ci+i(m,n)} 


(25) 


(26) 


m=0 n=0 


where Imf} represents the imaginary part of a complex quantity. The subsequent calculation of A* and Ay from e* 
and e is straightforward as indicated in Eqs. (23) and (24), except that the statistical averages are replacedby time 
averages performed by low pass filters. The calculated A* and A„ will be used to update the movement estimates 

through an accumulator, such that 

(27) 

(28) 


xi + 1 = xi + A* 
yi+i =yi + A„ 


The updated accumulator contents will be used to calculate the estimate <9 ro ,„.j+i and prepare the translated reference 
image for the next loop iteration at tj+i- 


5. CONCLUSION 

This paper describes an extended-image spatial acquisition and tracking scheme for an optical downlink, 
transfom-domain image correlation between the imaged downlink laser beam and the pnorly known £ 
intensity profile of a distant point source to determine the relative movement of the laser spot on the focal plane and 
the fine steering mirror to maintain accurate and stable optical pointing. 
based technique provides the optimal solution to spatial acquisition and tracking and m e^ected te perfom bet 
with low signal-to-noise ratios. The direct benefits from this fine-pointing technique mclude: (1) less burden on the 
telescope’s coarse tracking mechanism - a very cost-effective way to enhance pointing accuracy, and (2) possible 
simultaneous beam tracking and telemetry data detection on a single high-speed detector array minimizes th 
handover and possible alignment problems between acquisition and tracking operations. 






Preliminary numerical simulation using a computer-generated 2-D Gaussian intensity profile have shown that sub- 
pixel resolutions can be achieved by this correlation-type acquisition and tracking scheme in a high disturbance 
environment. However, more realistic simulation and demonstration to include small-sized focal plane array and real 
laser beam profile subject to various optical turbulence effects are required to determine the relationship between 
system performance and important specifications such as minimum loop update rate. 

There are several areas remain for further investigation. For example, the relative movement considered in this paper 
is limited to be rotation-invariant. The technique to compensate rotational movement encountered in image tracking 
is currently under study. The need for adaptation of reference image to compensate the line-of-sight optical turbulence 
is also an important issue to be resolved before successfully applying such a correlation-type technique. Furthermore, 
the assumption of independent white Gaussian disturbance in each pixel does not address many different types of 
disturbance. For example, sub-pixel scanning technique for resolution enhancement and various atmospheric effects 
can cause spatially correlated disturbances which will be an interesting area for future studies. 
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